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ABSTRACT 

Einstein's equations for stationary axisymmetric fields are reformulated as the equa- 
tions for affine geodesies in a two-dimensional space. The affine coUineations of this space 
are investigated and used to relate explicit solutions of Einstein's equations with different 
physical properties. Particularly, the solutions describing the exterior fields of a dyon and 
a slowly rotating body are discussed. 

PACS No. 04.20.-q, 04.20.Fy 



1. Introduction 

To simplify the structure of Einstein's equations, it is usual to postulate the existence 
of one or more Killing vector fields in the spacetime under consideration or, in less technical 
terms, the independence of certain coordinates. In a more general sense, the omission of 
the coordinates can be regarded as a special case of the Kaluza-Klein approach. Indeed, 
to investigate solutions with two Killing vectors in a systematic fashion, we can consider a 
Kaluza-Klein type reduction of Einstein's theory to two dimensions [1] . The dimensional 
reduction just amounts to dropping, for all the fields in the spacetime, the dependence on 
the coordinates that can be associated with the Killing vectors. 

In this work, we are concerned with a different type of dimensional reduction in which 
the number of fields - in our case, the metric coefficients - is reduced to the minimum 
necessary for describing the spacetime. This reduction occurs at the level of the Einstein- 
Hilbert Lagrangian and consists in dropping the terms that can be represented as total 
divergences, and rearranging the non-ignorable terms so that the Lagrangian becomes 
two-dimensional. The spacetime coordinates are absorbed into certain differential opera- 
tors that act on the remaining metric coefficients, i.e., the metric coefficients become the 
coordinates of 

the reduced space. This idea is in the spirit of the construction of the 
superspace studied by de Witt and others [2], where each point is a 
space. In our case, each geodesic defines a solution to the 
Einstein equations. 

Neugebauer and Kramer [3] introduced the abstract potential space 
determined by the Einstein-Hilbert Lagrangian of Einstein-Maxwell fields, and in- 
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vestigated the field equations which are derivable from a minimal surface problem in the 
potential space. In a recent work [4], we showed that canonical transformations can be 
used to reduce the dimensionality of the 

potential space, and Einstein's equations coupled to any matter field are 
equivalent to the geodesic equations in a two-dimensional space. 
As it is known, there exists a great deal of hypersymmetry in bidimensional physics, 
that is, supposedly unrelated problems happen to be the very same one, once a type of 
conformal transformation is given [5]. In our geodesic problem, 

this is equivalent to relate, or generate, new solutions to the affine 
geodesic motion, which, again refrased in the Einstein equations, translates into the 
generation of solutions. 

We will focus attention on stationary axisymmetric solutions to the Einstein equations. 
In Section 2, the corresponding field equations are derived from a two-dimensional metric 
Lagrangian and it is shown that they may be interpreted as the equations for an affine 
geodesic. We then investigate the equation for affine coUineations and present the general 
solution for the special case of a symmetry vector that depends on the coordinates only. 
Section 3 contains a solution generated by applying three different transformations on the 
Chazy-Curzon metric. We study the properties of this solution and show that it may be 
interpreted as describing the exterior field of a gravitational dyon. Section 4 is devoted to 
the derivation and study of a solution which contains the parameters necessary to describe 
the field of a slowly rotating mass. 

2. Field equations and affine coUineations 
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Consider the general stationary axisymmetric line element in Weyl canonical coordi- 
nates 

ds"^ = e^^^idt - ujd(l)f - e-2^[e2^(dp2 + dz"^) + p'^dcj)^] , (1) 

where ip, uj, and 7 are functions of p and z only. If uj =const., Eq.(l) leads to the special 
case of static axisymmetric fields. The calculation of the corresponding scalar curvature 
leads to the Lagrangian 

^ = ^(^p + ^l) + Mi^PP + i^zz - ipp - izz - V'p - '^D + sV'p , (2) 

which generates the usual Einstein equations. We proceed to give the main steps to 
construct from Eq.(2) another Lagrangian which generates the equation for affine geodesies. 
Introducing the differential operator D = (dp, dg), Eq. (2) can be written as 

L = 2DpD-f + {Diof - 2p{D'il;f . (3) 

2p 

In obtaining Eq.(3), we made the substitution pD'^B = D{pDB) — DpDB and neglected 
the total divergence terms. The Lagrangian (3) may now be interpreted as describing a 
kinematic system defined in the three-dimensional "space of metric 

coefficients" with generalized "coordinates" '0, a;, and 7 that depend on p and z. 
Consequently, p and z may be considered as quantities used to parametrize the coordinates. 
Equation (3) 

shows that in this special case the Lagrangian explicitly depends on the 
parameter p. Now, since 7 and uj are cyclic coordinates of 

the Lagrangian (3), it is convenient to use the Routhian R obtained from the La- 
grangian L by means of Legendre transformation acting on the cyclic 
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coordinates only, i.e., 



Here 11^^ is the canonically conjugate "momentum" associated with the generahzed coor- 
dinate U). 

Note that the conjugate momentum 11^ (as well as 7) does not enter the Routhian (4) 
at all. As a consecuence, it can be shown that the metric function 7 is determined by two 
first order partial differential equations that can be integrated by quadratures once and 
uj are known [6]. 

It follows from Eq.(4) that is a "constant of motion" (i.e., DH^ = 0) in space 
of metric coefficients. Using this fact, we can define an additional differential operator 
D = {—dz, dp) such that DD = 0. Introducing a function Cl by means of the relationship 

= p-^e^^Dw = Dfl , (5) 

the Routhian (4) becomes 

R=lpr'[{Dnr + {Dff] , (6) 

where / =exp(27/j). In this way we have obtained a new Lagrangian, Eq. (6), which 
corresponds to the squared line element of an abstract space described by / and Q and 
endowed with a two-dimensional metric tensor conformal to the Euclidean one, i.e.. 

The metric functions / and O plays the role of coordinates, and the operator D may be 
interpreted as the derivative with respect to the aflfine parameter (s) used to parametrize 
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the coordinates / and fl. Consequently, the spacetime coordinates turn out to be affine 
parameters on the space of metric coefficients. 

Therefore, Eq.(6) determines a metric Lagrangian that exphcitly depends on the affine 
parameter p. 

The geometric non-vanishing quantities associated with the metric (7) are, 
up to symmetries, 

ri — ri — r2 _ ^ 

^ 11 — ~^ 22 — ^ 12 — ~ J) 

-^212 = y2 ' (^) 

-Rii = -R22 = -p , 

and the scalar curvature is 4/p. 

The Euler-Lagrange motion equations obtained from the Routhian (6) are 

D^f - r\Df - + p-^DpDf = 0, 

(9) 

D'^n - 2f-^DfDQ + p-^DpDQ = 0, 
which are the same principal equations which follow from R^i, = 0. For completeness, we 
mention that taking E = f + iQ, the Routhian (6) can be rewritten as 

R= DE DE* , (10) 

where an asterisk represents complex conjugation. The variation of Eq. (10) with respect 
to E or E* leads to the Ernst equation [7] 

{ReE)AE= (DE)'^, with AE = D^E + p-^DpDE . (11) 

Since the starting Routhian (6) depends explicitly on the parameter p, Eqs.(9) coincide 
with the equations for an affine geodesic 
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D^X" + V\^DX^DX'' = X{p)DX'' . 

An affine geodesic is therefore a two-dimensional curve X^ — (/, fi) with tangent vector 
DX°^ satisfying Eqs.(9) for the real function A(p) = — 1/p. The existence and uniqueness 
of solutions of Eqs.(9) follow from theorems on systems of differential equations. Let Xq 
be a point of the space generated by gab and DXq the value of the tangent 

vector at this point. Then, there exists a unique, up to a change of the parameter, 
maximal affine geodesic X'' such that X'*(0) = Xg and DX''{0) = DXg (see, for in- 
stance, Ref. [8]). This result represents an alternative proof of the fact that a stationary 
axisymmetric solution is uniquely determined by its values on the axis of symmetry [9]. 

Let X^ and X2 represent two different affine geodesies. Our goal is to find transforma- 
tions that relate with and may be used to generate new solutions from known ones. 
The existence of this type of transformations cannot be assumed a priori and it depends 
on the symmetry properties of the underlying equations as well as on the explicit form of 
Xi and X2 . In the four-dimensional spacetime, transformations generating new solutions 
have been extensively studied and applied to diverse problems [10]; here, we first carry 
out the dimensional reduction and then reduce the problem to that of affine geodesies. To 
begin with the study of the transformations relating two different solutions of Eqs.(9), we 
consider the simplest case of a linear infinitesimal transformation. 

For the general affine geodesic equation, 

^2j^a ^ Y^^DX^DX" + g''^D{gi,c)DX'' = , (12) 
an infinitesimal transformation 

X"" ^ X"" = X"" + 6 7]"" , (13) 
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is said to be a symmetry transformation, that is, maps solutions into solutions, to order e, 
if the symmetry vector 77" satisfies the condition 

+ R%,dDX^DX'^V'' - {r\,),sDX^r^- = , (14) 
where D is the total derivative operator on shell 

^ = ^ + ^^"af^ - [^"bcDX'DX^ + 9'''Digu)DX^]-^ , (15) 

and s is a parameter along the geodesies. Moreover, R\^^ is the Riemann tensor of the 
space of metric coefiicients and period stands for 

partial derivative. If the metric Lagrangian is independent of the parameter s, Eq.(14) 
reduces to the equation of affine coUineations >Cr°'^^ = 0, where C is the Lie derivative along 
a curve with tangent vector ry", which is equivalent to the equation of the geodesic devi- 
ation for the connecting vector 77". Hence the geometrical basis of our approach becomes 
plausible. A family of solutions of Einstein's equations is equivalent to a congruence of 
geodesies in the space of metric coefficients. If r]°' is a vector connecting two 

neighboring geodesies at a given point, then the condition for rj'^ to remain a connecting 
vector at any other point of the space of metric coefficients, at which the 

geodesies are well-defined, is that it must satisfy the equation of geodesic deviation. 

If the symmetry vector r]"" is just a function of the parameter s and 

the coordinates, then the symmetry equation (14) can be rewritten as 

V%s + 2«J;6i^X^ + (r/ffe, + R%,a rj'')DX^ DX'= =0, (16) 

where a semicolon represents the covariant derivative associated with the metric gab given 
in Eq.(7). Notice that in the case that r]°' is just a function of the coordinates, even for a 
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metric depending on the non-affine parameter, the symmetry equation reduces to that of 
afSne colhneations. Consider this last case, ry" = 'r]"'{X^). Introducing the 
metric (7) into the symmetry equation (16), we get 

D^r]^ - 2f-\DfDr]^ - DVtDrf) + r]^f-^[{Dff - {D^f] + p-^DpDr}^ = , 

(17) 

D'^r]'^ - 2f-^{DfDrf + DQDr]^) + 2r]^f-'^DfDn + p-^DpDrf = . 
A detailed investigation of Eq.(17) shows that it possesses three independent solutions: 

Vi = (0, 1) , (18) 

V2 = if, ^) , (19) 

Vl = {f^,^^)- (20) 

Moreover, it can be shown that there are no affine eigencollineations, that is, the solutions 
(18-20) coincide with the Killing vectors of the metric (7). To find more general symmetry 
vectors of the potential space, it is necessary to consider the most general ansatz ry" = 
'q°'{s,X^,DX^). In this work, however, we want to focus attention on the symmetry 
vectors (18-20) and to show that even these simple vectors can be used to connect classes 
of solutions with different physical properties. 

We will now consider the type of solutions which can be generated by means of the 
vectors (18-20). Let ei, £2, and £3 be the parameters introduced by the symmetry vectors 
?7i, r]2i and ryf, respectively, according to Eq.(13). Acting on a seed solution (/, fi), the 
vector r]i leads to the new affine geodesic f = f and Q' = Q + ei. According to Eq.(5), this 
is equivalent to adding a constant coq to the metric function u. Obviously, this symmetry 
transformation is trivial since a coordinate transformation of the form t' — t — ujocp in the 
line element (1) absorbes the new term. Physically, this is equivalent to the introduction 
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of a rotating frame for the line element (1). Similarly, it is possible to show that the 
parameter e2 associated with the symmetry vector 773 can be absorbed by means of a 
rescaling of coordinates. The only non-trivial symmetry vector is 773 and it can be used 
to generate new solutions of the form 

/' = /(I + 6317) , n' = n+'-^{n'-f) . (21) 

Although, when acting alone, the symmetry vectors ry" and ryf are trivial, we will see 
below that they are helpful when used together with 773 to generate non-trivial solutions. 

Note, moreover, that the corresponding parameters ei and €2 can take any real value 
because they do not enter the symmetry equations at all. That is, 

putting the infinitesimal transformation (13) with rji and ryf > one sees that the result- 
ing equation is identically satisfied regardless of the values of the parameters ei and €2, 
respectively. 

3. Exterior field of a gravitational dyon 

The interest in monopole structures has rapidly increased during the past few years 
due to their discovery in generalizations of the standard model of particle physics. Magnetic 
monopoles were first introduced by Dirac [11] in electrodynamics to symmetrize Maxwell's 
equation in a direct way. Certainly, the most important consequence of the existence of 
magnetic monopoles is the quantization of electric charge. Most grand unified theories 
possess t'Hooft-Polyakov monopoles [12]. In general relativity there exist two different 
sorts of monopole structures: a magnetically charged black hole and a gravitational dyon. 
In fact, the magnetic black hole is the magnetic counterpart of the electrically charged 
black hole 
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described by the Reissner-Nordstrom metric, and is related to it by a duality rotation. 
A magnetic black hole can also be interpreted as a magnetic monopole with mass greater 
than a determined critical value [13]. 

A gravitational dyon is a hypothetical object the existence of which follows from the 
relativistic character of gravitation. In Newtonian theory, the only source of gravitation 
is the mass. In contrast, general relativity predicts that mass as well as rotation are sta- 
tionary sources of gravitational interaction. This leads to the well-known analogy between 
relativistic gravity and electromagnetism. The gravitational field generated by a distribu- 
tion of mass turns out to be analogous to the electric field, and the field of an angular 
momentum current presents characteristics similar to those of a pure magnetic field. For 
this reason, the field generated by an angular momentum current is called "gravitomag- 
netic" field. For this analogy to be complete, it is necessary to require the existence of a 
"gravitomagnetic monopole" as the counterpart of the magnetic Dirac monopole of electro- 
dynamics. A gravitational dyon is thus a mass endowed with a gravitomagnetic monopole. 
In this section, we will investigate solutions that can be generated from a static seed metric 
by means of a combination of symmetry transformations, and may be used to describe the 
exterior field of a gravitational dyon. 

To give a correct interpretation of the solutions presented here, we will use a 
coordinate-invariant method based upon the investigation of the relativistic multipole 
moments for asymptotically fiat solutions, according to the definition proposed by Geroch 
and Hansen [14] . We now proceed 

to derive the solution for a gravitational dyon. If we consider a static asymptotically 
fiat solution (/, f2 = 0) as seed metric and apply to it the symmetry transformation associ- 
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ated with the vector r/f , we obtain a stationary solution with f = f and 0,' = —€3/^/2. It 
can be shown that for any given asymptoticaUy flat / the new solution does not satisfy the 
condition of asymptotic flatness a la Geroch-Hansen [15]. Consequently, it is not possible 
to covariantly interpret the solutions generated by this type of transformation. To avoid 
this difiiculty, we use a combination of three different symmetry transformations (18-20). 
To the seed static solution / we flrst apply the symmetry vector rjf with parameter ei. 
The resulting solution is then used as seed solution for a transformation with the vector 
r]2 and parameter 62, and, flnally, we apply the symmetry vector 773. The new solution can 
be written as 

/'=(l + e2)/[l + eie3(l + e2)] , (22) 

and 

^^'=(l + e2)[ei-|(l + 62)(/^-e?)] . (23) 

It is now necessary to choose the parameters introduced by the symmetry transforma- 
tions such that the new solution becomes asymptotically flat. This condition leads to the 
relationships 

where 62 is a negative constant deflned in the interval €2 G (— 2, 0)\{ — 1}. As we mentioned 
at the end of section 2, the parameters ei and €2 do not need to be inflnitesimally small. 
Consequently, they can be chosen such that Eq.(24) is satisfled and €3 becomes inflnitesi- 
mally small as required by the transformation law (13). In fact, even for very large values 
of ei, €3 remains inflnitesimal and €2 remains in its domain of deflnition. 
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To analyze a concrete solution, we have to specify the asymptotically flat seed metric. 
Consider the Chazy-Curzon metric [16] 

/ = exp(-2m/r) , = + , (25) 

where m is a positive constant. The new solution is then given by Eqs.(22), (23) and 
(25). Choosing the new parameters according to Eq.(24), we calculate the corresponding 
Geroch-Hansen multipole moments and obtain 

Mq = m , Jq = —mes . (26) 

There are higher mass multipole moments M„ which corresponds to the axisymmetric 
mass distribution of the source, and higher moments for the angular momentum current 
Jn which, however, can be neglected since they are proportional to €§. Equation (26) 
shows that this solution represents the gravitational field of a body with mass m and 
gravitomagnetic monopole —me^. Hence, the new parameter €3 may be interpreted as 
the specific "gravitomagnetic" mass which may be positive as well as negative. The total 
"gravitoelectric" mass of the seed solution has not been affected by the action of symmetry 
transformations. For the sake of completeness, we present the metric functions of the new 
solution: 

/' = exp(-2m/r) , J = -2me3(l + e2fz/r , 7' = -m^p^jr'^ . (27) 

Finally, we would like to mention that using the Schwarzschild metric as starting 
solution, it is possible to generate the linearized Taub-NUT (Newman-Unti-Tamburino) 
solution which is also a candidate for describing the exterior field of a gravitational dyon. 
In general, it should be possible to find other solutions which, being different from the 
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Taub-NUT metric or the one presented here, present similar properties and hence might 
be used to describe a dyon. They all could differ only in the set of multipole moments 
higher than the monopole one; that is, there may exist different distributions of mass 
possessing the same gravitomagnetic monopole structure. 

4. Field of a slowly rotating mass 

For the study of the gravitational field of astrophysical bodies like stars and planets, it 
is necessary to investigate solutions which possess a set of mass multipole moments as well 
as a set of gravitomagnetic moments representing the rotation of the source. In contrast 
to the solution presented in the last section, a solution with realistic rotational properties 
may have only gravitomagnetic multipoles higher than or equal to the dipole one. In this 
section we derive a solution wich satisfies this condition. 

Consider any stationary seed solution (/, Q) satisfying the conditions of asymptotic 
flatness. As we have done in section 3, we apply three consecutive symmetry transforma- 
tions according to Eqs.(13) and (18-20). The new solution is then given by 



In general, this new solution is not asymptotically flat. However, if we demand that the 
parameters ei and €2 satisfy the relationships (24), asymptotic flatness is conserved and 
the resulting solution can be written as 



/' = (l + e2)/[l + e3(l + e2)(l^ + ei)] , 



(28) 



n'={l + e2) n + ei 



(29) 



(30) 



(31) 
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The calculation of new solutions does not present any difficulties. We will present here 
only one solution which illustrates our approach and can easily be interpreted. Consider 
the seed solution [17] 

^ -l + al{y^ - 1) ^ 2ai{x + y) 

^ {x + 1)2 + aliy - 1)2 ' {x + 1)2 + aliy - 1)2 ' ^ ^ 

with 

a;=^(r++r_), y=^(r+-r_), r| = + (2; ± m)^ , 

where m and ai are constants. To illustrate the effect of symmetry transformations, we first 
analyze the seed solution (32). An investigation of the corresponding multipoles show that 
there are gravitoelectric as well as gravitomagnetic monopole and dipole moments. Due 
to the presence of the gravitomagnetic monopole and gravitoelectric dipole, this solution 
cannot be considered as a candidate for the description of the gravitational field of any 
astrophysical object. Hence solution (32) is of no interest from a physical point of view. 
However, if we apply three different symmetry transformations to solution (32), its physical 
meaning can totally be changed. In fact, putting Eq.(32) into Eqs.(30) and (31), and 
calculating the relativistic multipole moments of the resulting solution, we see that all 
undesirable multipole moments vanish if cki is assumed to take the value 

ai = -€3(1 + €2^ . (33) 

Then, the only nonvanishing multipoles are 

Mo = m , and Ji = €3(1 + e2fm . (34) 

The last equation shows that the total mass of the body is given by m and that only 
the gravitomagnetic dipole moment survives in accordance with the dipole character of 
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rotation. The angular momentum per unit mass is given by €3(1 + €2)^ and can be positive 
as well as negative, corresponding to the two possible directions of rotation of the source 
with respect to the symmetry axis. Consequently, the new solution may be interpreted as 
describing the exterior field of a slowly rotating mass. Using Eqs.(30)~(33) and (5), the 
calculation of the metric components leads to 

/'=^, c.' = 2m63(l + 62f^, y=^ln4-^. (35) 
X + 1 X — 1 2 x^ — y 2 

This is equivalent to the Lense-Thirring metric [18], the physical meaning of which has 
been investigated by using other approaches and coincides with that we have obtained 
above by just analyzing the corresponding multipole moments. 
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5. Conclusions 

We have presented a different way to view to Einstein's 
equations, mainiy as geodesic motions in a space wliere tlie metric 
coefficients of tfie spacetime play tfie role of coordinates. The 
approach presented here for the axisymmetric stationary case can be 
generalized to any spacetime and even to the non-vacuum cases. We 
have reasons to believe that a dimensional reduction, via canonical 
transformations of the Hamiltonian, can always be made in this space 
of metric coefficients, and field potentials, such that the dynamical 
problem reduces to study the geodesic motion in a two dimensional 
manifold. 

As an application of this point of view we studied the symmetries of 
the geodesic motion for the space associated with axisymmetric 
stationary gravitational fields, and were able to generate some solutions, 
whose 

implications are currently under study. Nevertheless, we want to 
stress the fact that the idea presented here is not only a method for 
generating solutions but more than that, a different point of view to 
work with Einstein's equations. 
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